Origin of the Optical Gap in Half-Doped Manganites 



(N 
O 
O 

m 



Mario Cuoco,^ Canio Noce/ and Andrzej M. Oles^''^ 
^I.N.F.M. -Unitd di Salerno, Dipartimento di Fisica "E.R. Caianiello", 
Universttd di Salerno, Via S. Allende, I-84O8I Baronissi, Italy 
^'Marian Smoluchowski Institute of Physics, Jagellonian University, Reymonta 4, PL-30059 Krakow, Poland 
Max-Planck-Institut fiir Festkorperforschung, Heisenbergstrasse 1, D-70569 Stuttgart, Germany 

We have analyzed the coexisting charge and orbital ordering in half-doped manganites using a 
model which includes Coulomb and Jahn- Teller orbital polarization interactions. Most surprisingly, 
the gap in the optical conductivity is reduced by both on-site and intersite Coulomb interactions, but 
increases and explains the experimental results when the Jahn- Teller terms with orbital polarization 
are considered. The origin of this behavior is explained within a molecular model which arises in 
the limit of extreme topological frustration, when single electrons are confined to molecular units 
consisting of three orbitals. 
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I. INTRODUCTION 

The colossal magnetoresistance (CMR) manganites 
form a class of materials in which the problem of the 
interplay between the spin, orbital, and charge degrees 
of freedom is believed to lie at the heart of understand- 
ing their physical properties. ^'^ The parent compound, 
LaMnOa, is a Mott insulator due to large local Coulomb 
interaction U, with a single electron localized in degen- 
erate eg orbitals at Mn'^+ ions. The magnetic and or- 
bital ordering follow then from superexchange interac- 
tions which stabilize antiferromagnetic (AF) order. The 
doping by holes adds the charge degree of freedom in 
systems like Lai_a;Sra;Mn03, and triggers a transition to 
ferromagnetic (FM) state, in which the kinetic energy 
of the doped eg holes is optimized. This mechanism, 
called double exchange, was postulated five decades ago 
by Zener.'' Yet, the AF order is favored again at half dop- 
ing (x = 0.5),^^^ indicating that the AF superexchange 
and the FM double exchange frustrate each other, as in 
the FM metallic phase. ^ 

The half-doped manganites are very peculiar, and ex- 
hibit a so-called charge exchange (CE) AF ordering, 
which involves FM zigzag chains that are staggered anti- 
ferromagnetically; it has been observed experimentally 
by neutron scattering^'^ and by X-ray diffraction, 
both in the one-plane Lao.5Sri.5Mn04,^'^*' and in cubic 
Ndo.sSro.sMnOs manganites. ^'^^ The mechanism of sta- 
bility of the CE phase is still under debate. The hop- 
ping between degenerate Cg orbitals along zigzag chains 
shows a topological frustration^^ which leads to a band 
insulator. Although the CE phase is then favored, the 
large on-site Coulomb element U destabilizes it for realis- 
tic parameters. Therefore, one expects that either inter- 
site Coulomb interactions, ^^'^^ or the coupling to the lat- 
tice which leads to the Jahn- Teller (JT) distortions, ^^'^^ 
might play an important role. 

We argue that, similar as for the FM metallic phase, 
the optical experiments may provide the clue to the 
proper understanding of the CE phase within a truly 



relevant microscopic model. The experiments for the 
Lao.5Sri.5Mn04 compound^^'^^ show peculiar features by 
decreasing temperature from 300 K to 10 K. Above the 
charge /orbital transition Too — 217 K, the midinfrared 
spectrum is broad with an optical gap of the order of 0.2 
eV and a width of about 1.5 eV. It centers around the 
main peak at ~ 0.8 eV. Below Tco one observes a rigid 
shift toward high energy and a narrowing of the spec- 
trum. At the same time, the main peak shifts towards 
the bottom of the spectrum, giving rise to an asymmetric 
profile of the optical curve and an increase of the spec- 
tral weight. Finally, as the temperature is lowered below 
the Neel temperature Tat ~ 110 K, a further slight and 
smooth increase of the optical gap and narrowing of the 
main peak are observed. 

On the theoretical side, there have been only few re- 
ports on the optical properties of the charge and orbital 
ordering in half-doped manganites, ^^^^^ based on the elec- 
tronic structure calculations. In fact, these studies can 
serve as a starting point for the analysis of the corre- 
lated character of the optical excitations in the case of 
Cg electrons which are coupled either via the Coulomb 
interaction, or via the JT distortions. 

In this paper we introduce a new microscopic model 
for half-doped manganites. Next to the on-site Coulomb 
repulsion, it includes the JT term responsible for the eg 
orbital polarization at Mn^+ ions, induced by the JT 
distortions around the occupied sites when holes occupy 
the neighboring Mn^+ ions.^'^^ We study this model by 
exact diagonalization (ED) using Lanczos algorithm, 
and investigate: (i) the microscopic origin of the CE 
phase, its stability and properties in presence of the JT 
terms, (ii) the nature of the optical excitations as due to 
transitions between the correlated states of eg electrons 
at Mn(3d*) sites, and (iii) the reasons for the peculiar 
changes observed in the optical spectra under increasing 
temperature. ^^'^^ The unexpected result of the reduction 
of the charge gap due to electron correlations is explained 
as a consequence of the subtle role played by the topolog- 
ical frustration. This behavior turns out to be relevant 
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in explaining the pseudogap like features observed in the 
optical spectrum of the high temperature phase, where 
the distortions arc weak and the charge and orbital cor- 
relations are only at short range. Furthermore, we show 
that the large optical gap observed in Lao.5Sri.5Mn04 
at low tcmperatiire^^ cannot be simply explained within 
neither the electronic structure theory,^^ nor using a cor- 
related band insulator, but instead by including the JT 
interaction. 

The paper is organized as follows. In Sec. II we present 
the model Hamiltonian (Sec. II. A) and discuss the sta- 
bility of the CE phase, using as stabilizing mechanism 
either the JT term, or the intersite Coulomb interaction 
(Sec. II. B). Next we elucidate the frustration present in 
the kinetic energy due to the orbital phase factors, and 
introduce a so-called molecular model (MM) in which this 
frustration is maximized (Sec. II. C). The optical spectra 
are determined and analyzed both for the CE phase and 
the MM in Sec. III. A. The optical gap is extrapolated 
to the thermodynamic limit using finite size scaling (Sec. 
III.B). Finally, we discuss the obtained optical spectra in 
a broader context of the existing experiments, and sum- 
marize our results in Sec. IV. 



II. THE MODEL AND ITS GROUND STATE 



A. Importance of the JT coupling 

The model we adopt follows from the experimental ev- 
idence in manganites on the importance of the orbital 
degrees of freedom, the local Coulomb interactions oc U, 
and the JT distortions oc Ejt- The local Coulomb repul- 
sion U reduces the probability of double occupancies of 
eg orbitals, while the Hund's exchange coupling Jh be- 
tween the eg and the t2g electrons is responsible for the 
high-spin states, realized both at Mn"^+ and Mn'*+ ions. 
This term will not be explicitly included, as we assume 
that the eg and t2g spins are aligned, which corresponds 
to the frequently considered limit of Jh — > oo. 

We consider (spinless) Cg electrons moving along a FM 
zigzag chain of the CE phase shown in Fig. 1. Along each 
chain two kinds of nonequivalent sites alternate: bridge 
{i e B) and corner (j e C) sites (see Fig. 1). The 
polaronic model. 



Hpoi = Ht + Hu + HjT, 
consists of three terms, 



(1) 



Ht= J2 [(-l)^-ti('/')^'ls-.-i2(</>)6la,., + H.c. 



(2) 
(3) 



which stand for the kinetic energy (-ffj), the on-site (in- 
tcrorbital) Coulomb interaction U at corner sites {Hu), 
and the JT polaronic interaction (iJjx), respectively. 
The kinetic energy Ht describes the hopping oc f of Cg 



electrons between the directional 3x^ 



(3y^ 



or- 



bitals at bridge ions, and the x — y'^ {\x)) and 3z — r 
{\z)) orbitals at corner ions (see Fig. 1), with the corre- 
sponding electron creation operators: 6|, aj^, and a]^. 
The hopping elements are: ti{4>) = tshi(f), t2{4>) = 
tcoscf), and (/) = ^, with t standing for the hopping be- 
tween two identical directional orbitals along the bond 
(ij), i.e., 3.x^ — orbitals along the bond (ij) || a axis. 
The phase factors (— 1)^*^ follow from the alternating or- 
bital phases due to the zig-zag pattern, with Xij = 1 
or Xij = 0, where the hopping connects from the bridge 
position Bi or B2 with the corner site C, respectively. 




FIG. 1. Schematic picture of the charge and orbital (top 
panel), and of the spin configuration (bottom panel) for 
the CE phase. The arrows in the upper panel indicate the 
in-plane oxygon distortions around the manganese site at the 
bridge (61,62) and corner (C) position, which induce the cou- 
pling between the eg electrons and the JT disortions {qi}, as 
discussed in the text. The length of the arrows in the top 
panel is related to the amplitude of the distortion as due to 
the JT polaronic term. 



The JT orbital polarization interaction (4) splits off 
two eg orbitals at site «, if a neighboring site is occu- 
pied by a hole. This interaction has a similar form to that 
deduced for orbital polarons^® it favors the occupancy 
of the directional \ () orbital, oriented along the bond {ij), 
over the orthogonal to it planar |^) orbital (e.g., \3x'^ — r^) 
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over — z^) orbital for a bond along a axis). The classi- 
cal (phonon) variable qi stands for the respective oxygen 
distortions around an electron at site z, either a bridge or 
a corner position along the chain (Fig. 1), while g is the 
coupling constant between eg electrons and distortions of 
MnOg octahedra, and kjT is the spring constant for the 
JT mode distortions. It is convenient to introduce the 
JT energy, defined by Ejt = g'^ /kjT, which is naturally 
given by the scaling of qi = {g/kjT)qi, where g/kjx is 
the typical length scale for the JT distortion. By means 
of these substitutions, Hjt reads as: 



)(1 



(5) 



The JT term couples corner and bridge sites in an (a, 6) 
plane, either along the considered zigzag chain, or be- 
tween two adjacent chains. We minimized the total en- 
ergy (Hpoi) to eliminate the phonon variables {qi}, which 
gives the electronic Hamiltonian solved selfconsistently 
with: 



^{{ni(^-ni{){l-nj)). 



(6) 



The classical variables at bridge and corner positions 
are different, as suggested by different oxygen distortions 
around Mn ions.^^ If one considers the form of q-i, it is 
easy to recognize that the distortions are weak within this 
approximation, if the populations of the two orbitals are 
the same. Moreover, even in presence of a homogeneous 
orbital distribution, a small distortion pattern might be 
induced by the nearest neighbor orbital correlations. 



«rv y 



u \ 



FIG. 2. Schematic representation of the orbital structure 
in the MM. The chain is divided into molecular triads con- 
stituted of a bridge configuration (gray circle) together with 
an odd (circle with — ) and even (circle with -|-) superposition 
of the corner orbitals. The solid line shows the hopping pro- 
cesses inside the molecule, while the dashed lines indicate the 
intra- and inter-molecular interactions which result from the 
U and V coupling. 



different from that usually assumed, where the on-site Q2 
and mode are coupled to the local orbital density. -'^^ 
As regards the present model, we are assuming that a 
directional JT [Q^ like) mode, with its direction selected 
by the considered bond (ij), is activated by the pres- 
ence of a hole on the neighbor manganese, thus inducing 
a non-local correlated dynamics between the lattice and 
the charge degrees of freedom. As a direct consequence 
of this non-local coupling, one ends up with an effective 
directional Coulomb interaction between nearest neigh- 
bor manganese which strongly influences both the ground 
state and the properties of the low-energy excitations. 



B. Stability of the CE phase 

We begin with the investigation of stability of the CE 
phase with respect to the C phase, which consists of lin- 
ear FM chains along the h axis staggered in the a direc- 
tion, using the mean-field (MF) approximation. The C 
phase (01) is described by Hpoi (1) at (j) = |, with the 
hopping element t between directional — orbitals 
along each FM chain. While the on-site Coulomb element 
U increases the energy of the CE phase, it plays no role in 
the C phase, as the — orbitals are empty. Thus, the 
C phase is favored above [/ ~ 2.7t ([/ ~ 12.0t) in the MF 
(ED) method, if Ej^ = 0. We have verified that increas- 
ing Ejt induces in MF a charge disproportionation 5n 
in both phases, but this effect is more pronounced in the 
CE phase, lowering its energy especially when iJjx > t. 
Thus, we recognize that the JT term Tijx in Eq. (1) 
is the main stabilizing mechanism of the observed CE 
phase. The JT term induces a charge ordering (coexist- 
ing with orbital ordering) , which is in principle unlimited 
and may come close to 5n = 1, unlike the weak charge or- 
dering induced by the Coulomb interaction U, for which 
the order parameter cannot exceed 5n = 0.19.^'* Already 
at ?7 = the charge order is more pronounced than that 
induced by J/ —> 00 in a broad range of Ejt > 0.7t (Fig. 
3). At finite Ejt the charge disproportionation is only 
weakly enhanced by the on-site Coulomb element U . Al- 
though the form of the JT interaction term is different, 
our results agree qualitatively with the recent ab initio 
density- functional calculations.^^ 

TABLE I. Eigenstates and eigenenergies for a molecular 
unit i in the MM {(f) = ^) at U = Ejt = 0. The bonding, 
nonbonding, and antibonding states are labelled as \B}, \N), 
and 1^). 



state 


creation operator 




energy 


\B) 




-) 


Eb = -V2t 


\N} 






En ^0 


\A} 


4 = ;^&l-i(«I-i,+ + al+i, 


-) 


Ea = V2t 



We would like to point out that the form of the in- 
teraction between the JT modes and the eg electrons is 
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The CE phase ground state could also be repro- 
duced by intersite Coulomb interactions in the electronic 
model/^ 



K 



elec 



Ht + Htj + H^/ 



(7) 



with Hv = yJ2{ij) ntiij. Taking V ~ |-Ejt the models 
(1) and (7) lead indeed to similar MF phase diagrams. 

This might suggest that Tioicc could also provide a valid 
explanation for the observed CE phase. However, going 
beyond the MF one finds that the two models are fun- 
damentally different and lead to qualitatively different 
results, both for the ground state, and for the excita- 
tion spectra. 




1 2 3 1 2 3 
Ej-j/t Ej^t 

FIG. 3. Electron densities ric,: nt, Ux and Uz from toj} to 
bottom, as functions of Ejt, obtained using ED for tlie chain 
of L = 12 atoms at = (dashed lines) and a,t U = IQt 
(solid lines) in: MM (0 = f ), and CE phase (^ = f ). 



C. Molecular model 

Perhaps the most peculiar feature of the CE phase is 
the topological frustration of the kinetic energy which 
occurs due to conflicting orbital phases, and leads to a 
band insulator already at C/ = £^jt = 0.^® We empha- 
size that this frustration is maximized when the hopping 
elements with the opposite phases precisely compensate 
each other, i.e., ti — t2 — ^t, in the MM at 4> = j (see 

Fig. 2). Transforming then the orbital basis {a]^, a]^} to 
even and odd symmetry states at each corner site: 



,t 



V2 



(4 



(8) 



one finds that the hopping term simplifies to: 

Ht=tJ2 + fejai+i,- + H.c.) , (9) 



i=2n 



and represents a superposition of molecular units. Each 
of them consists of a central {i = 2n) bridge orbital b\ , an 
even symmetry + orbital at the preceding site i — 1, 

and an odd symmetry a\^^ _ orbital at the following site 
i + 1, with the eigenstates given in Table I. 
In the ground state at half doping. 



i^o>=n^iio)' 



(10) 



all Cg electrons arc confined to individual molecular units 
and occupy the bonding states, and the energy per site 
Emm = —Q.lQlt is significantly lower than in the CE 
phase, £^CE = — 0.695t. It is interesting to realize that 
although the electrons are now partly localized within the 
molecular units, the kinetic energy is at the minimum al- 
lowed in the subspacc of t\-\-t2 = ■ 

While the density distribution {ua} {a = b,x,z) is 
uniform for the noninteracting electrons, with the same 
density at bridge {rib = + n^) and corner (rig) posi- 
tions, increasing U or Ejt polarizes the system and in- 
duces a charge transfer towards bridge atoms. The effect 
of U is somewhat more pronounced in the MM than in 
the CE phase. This is best illustrated by considering the 
Gutzwiller projected state at U = 00,^'' which leads to 
the kinetic energy: 



(11) 



i=2n 



where a 



3:± 



(1 



rij zfi) arc the electron annihila- 



tion operators in the restricted space which does not 
contain double occupancies. By decoupling the quartic 
terms in Ht,^'^ bja^^^n^^^ ~ ^^S,±^"j,t) + (^Uj,±)^j,t' 



we have found a somewhat larger charge anisotropy Sn 
rib - (rix+nz) ^ 0.28 for the MM than 0.19 found before 
in the CE phase. A much stronger charge anisotropy, 
however, is induced in both cases by the JT term; with in- 
creasing Ejt it approaches gradually the localized limit, 
and at Ejt = 2t one finds already 5n > 0.65 (Fig. 3), 
with rix > riz due to the different hopping elements in the 
CE phase. The on-site Coulomb repulsion (3) is invariant 
under the transformation to {|-|-), |— )} orbitals and gives 
the interaction Uui^rij^ between two electrons from ad- 
jacent molecular units in the MM. It is quite remarkable 
that the effect of U in the ground state is weak in the 
regime of Ejt > t (Fig. 3). This demonstrates that the 
correlation effects are partly suppressed by the coupling 
to the lattice, similar as for the on-site JT term.^^ 



III. THE OPTICAL CONDUCTIVITY 

We now turn to the optical excitations and analyze 

the evolution of the optical gap A with the increasing 
Coulomb interactions and the JT coupling. The charge 
gap has been extracted by the calculation of the optical 
conductivity <t{lo) within a generalization of the Lanczos 
diagonalization technique. Since the system is in an in- 
sulating state, one can restrict the calculation of <j{ll)) to 
the finite frequency response given by the Kubo formula 



1 _ g-ZSi^ 
(j{uj) = Re 



f 

Jo 



rfre'--(j-(T)j), 



(12) 



where /3 = 1 /kgT and j stands for the contribution of the 
current operator along the FM zig-zag chain, obtained by 
the projection at the x{y) direction: 
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(13) 



Here Rim is the component of the vector which connects 
nearest neighbor sites I and m, tim is the hopping ele- 
ment between bridge and cornea' sites as defined above, 
and a = {x,z), respectively. By applying an electric field 
along the zig-zag chain, shown in Fig. 1, one accelerates 
charge along the chain, either along the x or along y axis, 
depending on the bond direction. In this way the current 
includes the charge-transfer processes from Bi or B2 to 
C sites, and vice versa. 



A. Interaction effects in the optical spectra 
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FIG. 4. Optical spectra cr{u!) as obtained in ED for the 
chain of L = 12 atoms in the polaronic model (1) with tp = 
and for: (a) U = Ejt = 0, ED result (full line) compared with 
an exact spectrum for CE phase (dashed line); (b) U = lOt, 
Ejt = and Ejt = 1.8t (dashed and full line). 



First, wc investigate the charge transport along the 
chain by analyzing the position of the lowest energy ex- 
citation in the optical spectra, which allows to extract 
the dependence of the optical gap on the value of the JT 
coupling constant and the electron-electron interaction 
parameters. Noninteracting electrons give a band struc- 
ture gap Ao = t in the CE phase (shown in Fig. 4),^"^ and 
Ao = \/2t in the MM. The on-site repulsion U typically 
enhances the energy of local excitations in the correlated 
systems, and increases the optical gap.^° Thus, it was 
unexpected that the gap first decreases linearly with U 
at E,^T = 0, and next saf,urates at finite A ~ 0.7Ao in 
the regime of t/ ^ t, both in the CE phase and in the 
MM (Fig. 5). This behavior can be understood, how- 
ever, by considering the lowest optical excitation in the 
MM for U o±t,\^i) = N|B^\^o), and expanding the in- 
teraction term oc Uni-nj+ using the Hubbard operators 
(X/s = A]Bj, etcetera), 

TV2{Xf'' + Xf ^ - Xf ^ - Xf «)] . (14) 

Indeed, the diagonal term oc X^^^ enhances the exci- 
tation energy by gJ/, suggesting that the gap A should 
increase. However, in spite of the electron localization 
within the molecular units, the above excitation is not 
localized, but can propagate along the chain due to the 
terms oc Xf ^X,^f in Eq. (14) [see Fig. 6(a)] with dis- 
persion — \U cosfc, resulting in a decrease of the optical 
gap calculated at band minimum, A ~ Aq — gf/. 
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FIG. 5. Optical gaps A found in ED for the MM (left) 

and the CE phase (right) as functions of electron-electron in- 
teractions (top): - U (circles), V (squares), and U -\- 2V for 
V = U/4, (triangles), and of Ejt (bottom) for: U = Q (cir- 
cles) and U = IQt (squaxes). Dashed line in top left panel 
shows the weak couphng expansion [Eq. (14)]. 

A similar weak coupling expansion performed for the 
intersite Coulomb term oc F in Eq. (7) gives rise in low- 
est order to the same optical excited state but it 
cannot propagate as now Heiec does not contain terms of 
the form oc X^^^X^^f . Therefore, acting on such states 
by Tioicc leads only to higher energy local excitations [Fig. 
6(b)], and explains why the optical gap increases first in 
the range of 2V < 5t, before it starts to decrease at higher 
values of V (Fig. 5). However, in the realistic regime of 
parameters (V -^U and V (xU) the gap would decrease, 
as dominated by the on-site term U. The gaps in the CE 
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phase follow the same trends when the electron interac- 
tions are increased (Fig. 5), with the optical gap being 
reduced to A ~ ().7t at U = lOt (Fig. 4). Thus, the weak 
propagation along the zigzag chain at </> = | does not 
change the generic picture of the gap decreasing with U, 
which emerges in the MM. 

In contrast, the interactions with the lattice described 
by the JT term lead to qualitatively different physics. 
Increasing Ejt decreases local potentials at bridge posi- 
tions in the MM, and thus the splitting between the \B) 
and \N) state and the gap value A increase. Although 
the gaps in the MM and CE phase are different in the 
realistic regime of Ejt ~ t, they approach the same limit 
A ~ 2EjT at large Ejt » t. Furthermore, even if the 
JT term contains an effective long-range Coulomb term 
mediated by the JT distortions, this contribution does 
not counterbalance the increase of the charge gap due to 
the splitting between the bridge and corner states. 
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FIG. 6. Schematic representation of the low-energy optical 
excitations induced in the MM by: (a) on-site U, and (b) 
intcrsite V Coulomb interaction. The lowcst-cncrgy excita- 
tion \4>i) is made starting from the ground state \4>o)- The 
terms oc U permute the bonding \B) and non-bonding \N) 
states on neighboring molecules and allow a propagation of 
the low-energy excitation along the chain (a), while the terms 
oc V provide only the conpHng of the low-cncrgy optical ox- 
citations to higher energy configurations which contain more 
than one nonbonding state in the chain (b). 

We also investigated the changes in the shape of the 
optical spectra <j{u)) with increasing interactions. First 
of all, the asymmetric spectrum for noninteracting elec- 
trons, with a large spectral weight at the low-energy fea- 
ture, is well reproduced (Fig. 4). Both increasing U and 
EjT lead to more pronounced asymmetry of the spectra. 
By analyzing the spectral weight distribution, we estab- 
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lished that electrons are excited from bridge (Mn 
corner (Mn^+) sites, and the x (z) character of the ex- 
citations is found predominantly close to the low (high) 
energy edge. The high energy part looses its spectral 
weight faster, and the distance between these two fea- 
tures decreases with increasing Ejt or U, resulting in 
much narrower optical spectra [e.g., the distance between 
two local excitations described by Eq. (11) is only 0.87i 
at [/ = oo compared with lAlt at U = 0]. 



B. Finite size scaling 



In order to extract information about the electronic 

system in the thermodynamic limit, we have performed 
a scaling analysis of the charge gap obtained for different 
values of the chain lengths (from L = 4 to L = 16 sites), 
and two representative sets of parameters with Ejt = 0: 
U = 0, 1, 3, and lOt, taking either V = ot V = 2t. 
We begin with the analysis of the influence of increasing 
U at V = [Fig. 7(a)]. In this case the charge gap is 
almost independent on the cluster size, which indicates 
that on-site excitations discussed in the previous Section 
indeed reduce the optical gap due to correlation effects. 
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FIG. 7. Optical gap A in the CE phase (in the units of 
the gap for noninteracting electrons Ao), found by means of 
Lanczos diagoualization for different chain size L = 4, • • • , 16, 
for V = 2t and: U = 0, 1, 3, Wt. The dashed hne in (b) indi- 
cates the linear interpolation for the critical size above which 
the charge gap saturates. 

To understand the results at V = 2t, let consider first 

the case oi U = 0. As shown in the Fig. 7(b), the behav- 
ior of the charge gap is monotonic up to L = 8 and then 
it saturates. One can observe that the only effect of V 
is to confine the charge in the ground state as the lowest 
charge excitation has a gap in the spectrum, and in addi- 
tion to localize the lowest optical excitation on a length 
which is smaller than L = 8, as the charge gap does not 
change above this critical size. In other words, the re- 
sults suggest that the quantum coherence of the lowest 
optical excitation extends only over a limited distance, 
so that when the cluster size exceeds L = 8 length, the 
increasing further cluster size does not affect the value of 
the charge gap anymore. 

Consider now the effect of the on-site interaction oc U. 
As one can see in the Fig. 7, increasing U reduces the am- 
plitude of the charge gap, and modifies the critical size, 
above which the gap saturates. Due to the limitation in 
the numerical computation, we can follow the evolution 
of the critical size (dashed line) only for U ^ t. Indeed, 
for large values ofU ^ t, one observes only a monotonic 
growth of A as a function of 1/L, without any sign of the 
gap saturation. We have extracted the behavior of the 
critical size vs U by means of a linear fit. One can argue 
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that the effect of U is to increase the locahzation length 
in the local optical excitation spectrum (LOES), which 
yields a larger value of chain size to get a saturation in 
the charge gap. This itinerant character induced by U is 
clear in the MM, where we have seen that finite U allows 
for a propagation of the LOES. 

By summarizing the results of the scaling, we deduced 
that the effect of V is to harden the charge fluctuations 
in the ground state and to induce a localization of the 
LOES on a length scale of about L = 8. On the con- 
trary, the U term induces an itinerancy in the LOES 
which manifests itself by an increase of the localization 
length in the LOES, and consequently of the critical size 
of saturation of the charge gap. The competition of these 
two counteracting effects turns out to give the reduction 
of the charge gap in the bulk limit only for a very weak 
intersite Coulomb repulsion, V <SCU. 

A further important conclusion from the presented 
analysis is that one finds a large range of variation of 
A/Ao. when one moves from a chain of L = 4 sites to the 
extrapolated bulk limit. This analysis might be relevant 
for real materials if one interprets the disordered high 
temperature phase of the half-doped manganites as made 
by a liquid of zig-zag FM chains with different length, so 
that the optical response has to be built by making a su- 
perposition of the response of chains with any size. We 
argue that the low-energy part of the spectra contains 
this fluctuating superposition of several chain contribu- 
tions of different size, and these short-chain contributions 
are suppressed when one enters into the ordered phase 
below Tco- 



IV. DISCUSSION AND SUMMARY 

We argue that the present model agrees qualita- 
tively with available experimental information, and our 
flndings reproduce the experimental observations in 
Lao.5Sri.5Mn04^^ and Ndo.5Ca1.5MnO3.31 In the high 
temperature regime when the JT distortions are weak, 
the optical gap is small and dominated by the effect 
of electronic correlations. When the temperature is de- 
creased below Tco I the spectral weight grows and shifts 
to higher energy, the optical spectrum is narrowed, and 
the asymmetry in the spectrum increases. All these fea- 
tures can he explained within the present polaronic model 
(1) as an indication of the onset of the orbital and charge 
ordering from the disordered state, described by finite 
values of local variables and thus increasing orbital 

polarization due to the JT terms oc EjTQi- Further shift 
of the peak feature up to A ~ 1.3 eV would follow from 
the increased {qi} below the magnetic transition, when 
the magnetic fluctuations are suppressed. The observed 
change of the gap A is reproduced by Ejt ~ 1.2t, which 
corresponds with t = 0.4 eV^ to the expected orbital 
splitting of ~ 0.3 eV.^^ 

It is worth pointing out that previous studies have 



discussed the optical conductivity in the CE phase as 
due to dipole transitions from Mn(eg) to Mn(4p) states 
by means of tight-binding approaches and first princi- 
ple band-structure calculations in the local spin-density 
approximation. ^3 The analysis of Ref. 23 included the ef- 
fect of the JT distortion as an effective energy level split- 
ting between the bridge and corner states. The proposed 
interpretation that the changes of the optical spectra are 
mainly due to the JT distortion on the Mn(eg)-Mn(4p) 
dipole transitions does not include the effect of correla- 
tions, and thus it fails to reproduce the experimentally 
observed modification of the spectral weight when the 
temperature is lowered. 

Our interpretation of the optical spectra is different. 
We have shown that the low-energy optical excitations 
are mainly due to the interband transitions between cor- 
related Mn(eg) states. The effect of the JT coupling 
cx EjT may be summarized as follows: (i) the uncorre- 
lated spectrum shifts toward higher energy, (ii) the spec- 
tral weight increases in the low-energy part, and (iii) the 
asymmetry of the profile grows as the spectral weight in 
the high energy shoulder gets suppressed. Particularly 
the second feature is strongly dependent on the form of 
the JT coupling. Indeed, by a simple level splitting of 
the bridge and corner states, one would just obtain a 
rigid shift of the spectrum, while the effective directional 
Coulomb interaction present in Hjt (5) is responsible of 
the increase of spectral weight in the low energy side. 

Finally, our results show that the shift of the main peak 
and the increase of its spectral weight come essentially as 
a consequence of the orbital polarization present in the 
JT coupling term. It influences the optical transitions 
between correlated Cg states. The pseudogap activation 
in the high temperature phase is mainly due to the ef- 
fect of Coulomb correlation which reduces the band gap. 
Moreover, as mentioned above, the low energy tail in the 
spectra observed in the high temperature regime, might 
originate from the superposition of the optical response 
of FM zig-zag chains with different length (see Fig. 7). 

In summary, we have shown that the subtle inter- 
play between the topology, the Coulomb correlations 
and the orbital polarization JT terms are responsi- 
ble for the temperature dependent optical properties of 
Lao.5Sri.5Mn04. Thus, we argue that these interactions 
play a crucial role for the observed coexisting orbital and 
charge ordering in half-doped manganites, which cannot 
be explained within purely electronic models. When the 
orbital polarization interaction is taken into account, the 
charge disproportionation turns out to be much stronger 
than that induced by U alone, but not as pronounced as 
in the original Goodenough picture. 
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